Giant flexoelectric effect in two-dimensional boron-nitride sheets 
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We find, with tlie use of first-principles calculations, that a single-atom-thick boron-nitride (BN) 
sheet exhibits an unusual nonlinear electromechanical effect: it becomes macroscopically polarized 
when bent out-of-plane. The direction of the induced polarization is in the plane of the film and it 
depends non-analitically on the corrugation wave vector k. The magnitude of the polarization can 
reach very high values in spite of being at least quadratic in atomic displacements due to BN sheets 
being able to tolerate large mechanical strains. The discovered effect can find many applications, in 
particular, in a new type of efficient and reliable nanogenerators. 

PACS numbers: 61.46.Np, 77.55. -l-f, 77.80.-e, 77.84.-s 



Three-dimensional (3D) bulk crystals can generate a 
voltage either in response to a mechanical strain diUi 
(piezoelectric effect [U [2]) or to a strain gradient didjui 
(flexoelectric effect [21 E]), where u{r) is the displace- 
ment vector, with I its Cartesian index, di the gradient 
operator. The flexoelectric effect is usually small and 
evades experimental detection unless large strain gradi- 
ents are externally imposed or artificially designed inho- 
mogeneous metamaterials are used [Sill])- It is commonly 
assumed that the effects nonlinear in diUi and didjUi on 
a polarization are negligible in bulk dielectrics. 

Contrary to 3D systems, sp^-bonded 2D crystals, like 
graphene and boron nitride (BN) [5l |5] , are able to sus- 
tain much larger structural distortions and, thus, exhibit 
new forms of electromechanical coupling. The BN sheet, 
for example, becomes pyroelectric when it is wrapped 
into a chiral or zigzag nanotube with the macroscopic 
polarization inversely proportional to the inverse square 
of the tube radius, l/R^, and directed along the tube 
[3 HJ E]. Formally, this effect can be considered as a 
quadratic flexoelectric effect since 1 /i?^ is only a partic- 
ular form of (didjUk)^. 

Here, we predict another unusual nonlinear electrome- 
chanical effect in 2D BN monolayer that has not been 
noticed so far: generation of a macroscopic in-plane po- 
larization in response to out-of-plane atomic displace- 
ments like the corrugation Uz(r) — ylsin(A; • r + ip), 
where k is the undulation wave vector, r is the in- plane 
vector, z— axis is perpendicular to the plane, and is 
some phase. Although such displacements have no in- 
plane components and produce zero net curvature, they 
nevertheless induce an in-plane macroscopic polarization 
P (k) , which is purely electronic in origin and related to 
the shifts of tt and a chemical bonds. It is remarkable 
that the induced polarization practically does not depend 
on the phase ip but strongly depends on wave vector k, 
and there are some ( "optimal" ) directions along which 
the absolute value |P(fc)| goes through a maximum. We 
show further that being at least quadratic in amplitude 
of atomic displacements A, the corrugation-induced po- 
larization is decomposed basically into a sum of two con- 
tributions associated with the quadratic terms {diUu) 



and (didjUf;)^, respectively. Whereas the flrst ("piezo- 
electric") contribution can be understood starting from 
the piezoresponse of an isolated planar sheet, the sec- 
ond ( "flexoelectric") term has the same nature as the 
curvature-induced polarization in BN nanotubes, where 

It should be stressed that the discovered effect bears 
no relation to the formation of the normal polarization 
due to bending of 2D systems like graphene flat sheets 
where P cx l/R |10[ lllj . From a more general point 
of view, the effect can be considered as a mechanically- 
induced improper ferroelectric phase transition, with A 
playing a role of an order parameter. The externally im- 
posed acoustic-phonon-like displacements Uz{r) reduces 
the initially nonpolar symmetry group to some lower po- 
lar subgroup, so that P (fc) is expanded in powers of A 
beginning with the term cx A^. Although P (k) is at 
least quadratic in A, the generated polarization can be 
very large provided that A is a fraction of the wavelength 
27r/fc. The effect may flnd a variety of practical applica- 
tions and, in particularly, for conversion of the ambient 
wave-like micromovements into electricity: this tantaliz- 
ing possibility is discussed below. 

Among the truly 2D crystals that have been recently 
obtained using the so-called cleavage technique [12] . 
boron-nitride monolayer is the only piezoelectric with 
wide band gap. Being partially ionic and partially cova- 
lent with sp^-bonding, a flat BN sheet has a remarkable 
mechanical flexibility associated with the ease of forming 
intermediate sp^^" bonds. Due to its D^h point symme- 
try group, it does not have a ground state polarization, 
although it exhibits piezoresponse with the piezoelectric 
tensor obeying the symmetry relations for the following 
non-vanishing components: 

^x,xx — ^x,yy — ^y,xy ^y.yx: (1) 

where x is directed along the symmetry axis of 2-th order 

US]. 

To describe periodically and commensurately distorted 
hexagonal BN sheets, we introduce the corrugation wave 
vectors k as 27re/A, where the unit vector e and wave- 
length A are expressed via the lattice vectors Oi and 



FIG. 1: Honeycomb BN lattice with basis vectors ai and 02. 
The X axis is chosen to be parallel to (01+02), whereas the 
y axis perpendiculat to it. The angle between the (01+02) 
and the chiral vector A is the polar angle 9 used in the paper 
to describe the the corrugation wave vectors k=2n / \{n,m). 
The example corresponds to the case of n = 3, m = 1. 



a2 of a 2D sheet in the following way: e = A/A, with 
X{n,m) — nai + ma2, A ~ a^/ n? + nm + , where n 
and m are integers, ai = 02 = a (Fig. 1). Note that the 
vector A(n, m) coincides with the chiral or circumferen- 
tial vector defined in the theory of carbon nanotubes [5]. 
The corrugation with the wave-length A(n, m) leads to 
a rectangular supercell whose translational vectors are 
A(n, m) and some translation vector T = Nai + Ma2, 
which is perpendicular to A. The vector T with length 
T = aV^/V2"+^VA/TM2 is completely defined by the 
vector A HH]: 

r VSX/p, N-M^ 3pq, , . 

\ VSX/ (3p) , N-M^ 3pq, ^ ' 



FIG. 2: A scheme showing the angular dependency of the 
polarization (depicted as arrows) induced by sinusoidal cor- 
rugations at a fixed " amplitude-to- wave-length" ratio Ak. 
While the fc- vector scans the angle 6 from to 2?!, the Carte- 
sian components of the polarization change approximately as 
(— cos 20, sin 29) , where 9 is the angle defined in Fig. 1. Green 
solid and blue dashed lines indicate the armchair and zigzag 
directions along which the polarization is pure longitudinal 
and transverse, respectively. 

process of corrugation: P — P^")- P^"^, with 

^'"'-^^E/./'K-i^l^-S); (3) 

where u^^^ is the periodic part of the Bloch functions 
and e is the electron charge. Note that the integral in 
this expression is taken over 2D Brillouin zone, so that 
the value of polarization has a dimension of a charge per 
unit length (electron charge per Bohr radius, e/ag, if 
atomic units are used). Alternatively, in the simplest 
case of isolated bands, the formula (l3| can be cast as 



where p is the greatest common divisor of N and M, 
and q is an integer. It is also convenient to define the 
vectors fc in the polar coordinate system, (fc, 9), where 6 
is the angle between the fc and (ai + 02) [IH]- We choose 
(ai + 02) in such a way that the shortest interatomic 
distance along this vector is passed from B to N atom 
(not from N to B, Fig. 1). Borrowing the terminology 
from the theory of nanotubes, we shall call the sine-wave 
distortions with 9 = zL ttI/3 the "armchair" -like and 
with 9 — 7r/6±7rZ/3 the "zigzag" -like, where I is an 
integer. All other distortions falling between the two will 
be called "chiral" . 

To calculate the corrugation-induced polarization, we 
use the Berry-phase approach [14]. Namely, we treat the 
electronic polarization as a geometrical phase accumu- 
lated by the occupied electrons in the slow (adiabatic) 
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where r„ is the center of the Wannier function (WF) 
corresponding to the occupied band n and S"^"-* — 
A T the surface area. Since the ionic contribution 
(e/S"^"^) Zii?-"'' to the total in-plane polarization 
does not change under the corrugations of interest, it 
can be completely neglected. 

The integrals P*^"' (3) were computed using ABINIT 
code [TS] with a 12x4x1 Monkhorst-Pack k-point grid, 
where the largest number 12 corresponds to the direc- 
tion of undulations A and 4 to the vector T, perpen- 
dicular to A. The flat and corrugated BN sheets were 
simulated by a slab-supercell approach with the inter- 
planar distances of 20aB and larger to ensure negligible 
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wave function overlap between the replica sheets. For the 
plane- wave expansion of the valence and conduction band 
wave- functions, a cutoff energy was used in the range of 
80-100 Ry, depending on the corrugation wave length A. 
The WFs and their geometrical centers r^"'' [see (4)] have 
been calculated with the PWSCF package [16] in three 
steps. First, we performed a self-consistent ground-state 
calculation and then a non-self-consistent calculation at 
fixed k-points for all the occupied bands. And finally, 
the obtained Bloch functions were transformed into max- 
imally locahzed WFs (MLWFs) as described in Ref. [Tfj . 
We used tt and cr orbitals as initial guess projection func- 
tions allowing to obtain converged result for MLWFs . 
There are Af 7r-derived and 3 cr-derived WFs, where 
Af is the number of B-N pairs per unit cell AxT. 

Before describing the results of ab-initio modeling, 
we shall derive a phenomenological expression for P{k) 
considering the BN sheet as a continuous elastic mem- 
brane, which appears to be very accurate. Obviously, 
in the long-wave length limit (fc — > 0) the induced po- 
larization can be expressed via 2D piezoelectric tensor 
components Sijk and pure shear deformations, because 
isotropic distortions do not contribute to the polariza- 
tion of a flat BN. Substituting the displacement field 
Uz (i") = j4sin(A; • r + ip) into the standard deformation 
tensor of a membrane 120] 



(5) 



1 / dui 
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where Xi — {x,y), one can easily find the net shear com- 
ponents 77 = (Ej-j. - eyy)/2 and 7 = e^,^ = as 



£l I cos 20 
Y \ sin 261 



(6) 



where ey = A'^k'^ /A is the net film stretch along the k 
associated with the third term in Eq. ([s]) and 9 is is the 
angle related to k the way shown in Fig. 1. The Carte- 
sian components of polarization in the long wavelength 
approximation become Px = 2ari, Py — — 2^7, where a 
is the 'clamped-ion' piezoelectric constant Cx.xx of a flat 
BN sheet. Using these components and relation (|6|, the 
vector of polarization can be presented as 



P{k) = 2e||(e|| cos 36* - e_Lsin36l), 



(7) 



where ey and ej^ are the unit vectors parallel and per- 
pendicular to k (e_L are defined so that z ■ [ej^ ^ 6||] > 0; 
z\\ [x X y\). It is clear that the vector P{k) simply ro- 
tates with 9 for the constant ey . 

Being accurate in the long- wave length limit (fc — > 0), 
the formula Q must be corrected for shorter wavelengths 
by supplementing the terms oc fc**. The first such term 
comes from the fact that the longitudinal component of 
polarization is tangent to the surface and tilted down 
and up relative to the initial {x,y)- plane. It, there- 
fore, should be projected onto the plane, which has 
not been taken into account in Eq. ([t]). To obtain 



this correction, one should average over A the quan- 
tity ^aA^/c^ey cos36'cos^(fe • r + ip) [cos'i9(r) — 1], where 
i9(r) is the tilting angle. Noticing that cosi?(r) — 1 w 
— ^A^k^ cos^{k ■ r + if) one can easily find the average as 
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yA'^k* en cos 39. 



The second correction accounts for the flexoelectric ef- 
fect, i.e. the appearance of the macroscopic polarization 
due to finite curvature like in BN nanotubes. By sym- 
metry, such a polarization is maximal in amplitude for 
the zigzag and minimal (zero) for the armchair direc- 
tions; this angular dependence can be simply described 
by sin 39. Besides, the polarization should be parallel to 
ej^, in close analogy with the tubes where the macro- 
scopic polarization is always along the tube. And, fi- 
nally, it should be proportional to the net inverse ra- 
dius of curvature squared, 1/i?^, instead of 1/i?^ in the 
tubes. Defining the principal local curvature along the fe 
as V^U2(r) = Ak'^ sin(fc • r + (p) ~ 1/R{r), one can easily 
find that TJW = = A^k*/2. Thus, the curvature- 
induced correction can be represented as j3 R~^e± sin 39, 
where /3 is the flexoelectric constant. 

By adding the above corrections to the expression ([7| , 
we come to the final result: 



P(fe) 



aeu 



-l3R-'^ej_ sin 361. 



(8) 



This formula is invariant with respect to rotation 9 
9 + 2tt/3, as one would have expected from the D^^ 
symmetry. The first term comes from the increase in 
surface area that accompanies the sine-wave atomic dis- 
placements, it represents the piezoelectric effect. The 
second term, on the contrary, is non-zero only due to fi- 
nite averaged curvature squared, obviously, it renders the 
flexoelectric effect. Since a < and /3 > 0, the extrema 
of the polarization are: 



Pniin = lakll 1 - :r£ 



zigzag, (9) 
armchair. (10) 



For all the other (chiral) directions, the amplitude falls 
within the above limits. 

The calculations show that for relatively small A/\ < 
0.2 the dependence P{k) does follow the continuous ex- 
pression (|8| very well. The mutual orientation of the P 
and k constantly evolves as the latter rotates through 
the angle 9 (Fig. 2]). P becomes pure longitudinal along 
the armchair fc— directions and pure transversal along the 
zigzag directions. Along the chiral fc-directions, it has 
both components. Being practically independent of the 
initial phase ip, the magnitude of polarization P slightly 
depends on fc (within 30%) and reaches its minima along 
the armchair and maxima along the zigzag directions. 

We have found a and /3 numerically to be — O.llSe/os 
and O.OAAeag, respectively. To check that the flexoelec- 
tric constant f3 is indeed related to the properties of BN 
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nanotubes, we considered a set of 'ideal' zigzag (n, 0) 
nanotubes cut out of BN sheets in such a way that the 
B and N ions fall on the cylinder with radius R, 2ttR co- 
inciding with the 'wave-length' A(n,0). These tubes are 
longitudinally polarized and their polarization is propor- 
tional to 1 /i?^ to a leading order in 1/R. The coefficient 
of proportionality, as we found from additional calcula- 
tions, is 0.042eaB, which is very close to the parameter /3 
(0.044eas). This means that the curvature-induced po- 
larizations in the corrugated BN sheets and in the BN 
nanotubes indeed have the same origin. 

To analyze the relationship between electronic struc- 
ture and P (fc), we notice that the WF in a flat BN sheet 
have roughly the character of the a- and 7r-bond orbitals 
(tJ and TT WFs). The tt— like WFs are centered exactly 
on the N atoms, whereas their a counterparts are some- 
where in the middle of the B-N bonds closer to the N 
atoms (Fig. 3C). The contributions of the tt and a WFs 
to the 'clamped-ion' polarization induced by an in-plane 
uniaxial stretching have the same sign, with the vr WFs 
dominating and amounting to as much as 85% of the 
total (see also Ref. [9 ). 

Although the corrugation of a flat BN sheet leads to 
the rehybridization effects like mixing of vr and tt* elec- 
tronic bands due to breaking the z — > — z mirror sym- 
metry it is still possible to classify all the WFs into 
TT and a types. Figs. 3 and 4. It might naively appear 
that in the corrugated structures the main contribution 
to -P(fe) would also come from the tt WFs. In fact, this 
is not the case. The contribution of a WFs generally 
increases, becoming strongly dominant in some special 
cases in corrugated films. 

Consider first the case of zigzag corrugation with n — 
4, m = -4, A/X = 0.1 and V5 = 0, Fig. 3. Here, the 
shifts of TT WFs are still dominant but not so apparently 
and iP-n/Pal is only about 2.5. In contrast to the flat 
case, the tt and a WFs have the opposite sign contribu- 
tions to P{k) exhibiting inhomogeneous shifts within the 
two-dimensional supercell A x T. Upon the corrugation, 
the TT WF centers shift relative to the initial positions as 
shown in Fig. 3D. Although the sum of the shifts along 
the y axis cancels within the supercell, the similar sum 
along the x axis is nonzero, which means the appear- 
ance of a macroscopic transverse polarization, P ^- k. 
Interestingly, the magnitude of the shifts is larger in flat- 
ter and more stretched areas, and the shifts themselves 
(along y) are always directed towards more curved and 
less stretched areas. Because of the charge transfer from 
the flatter to curved regions, the conduction band maxi- 
mum increases and the valence minimum increases in the 
curved areas, thus leading to a space modulation of the 
band gap with the wavelength A. 

Now, we turn our attention to another case, repre- 
sented by the armchair undulation with rt = 4, m = 4, 
A/\ = 0.1, and f = Q, Fig. 4. Here again, the tt and a 
WFs give the compensating contributions to the polar- 
ization. But now, in contrast to the zigzag case, a WFs 
dominate by a wide margin: \P^/Pt^\ ~ 4.5. All the WFs 



shift strictly either along or opposite the fe-vector (Fig. 
4C), so that the induced polarization is purely longitu- 
dinal. The dominant a WFs move towards the N-atoms, 
thus further increasing the polarity of the a B-N bonds. 
The magnitude of the shifts, as in the previous case, is 
larger in more flattened and stretched areas, but this 
works only in the second and fourth quarters of the A. In 
the first and third quarters, the shifts of a WFs are mod- 
est even in the well stretched areas. The charge transfer, 
therefore, increases the energy gap mainly in the second 
and fourth quarters at the expense of the first and third 
quarters of the wavelength. 

Taking a reasonable amplitude-to-wavelength ratio of 
A/\ ~ 0.2, it is easy to estimate from Eq. ([s]) that the 
induced polarization can be on the order of 4x 10~^e/oB. 
This quantity cannot be compared directly to that of 3D 
bulk materials because the latter has different dimension, 
e/a^. The comparison can be performed, nevertheless, 
if both 2D and 3D polarizations are expressed as the to- 
tal dipole per stoichiometric unit [3]; then the polariza- 
tion estimated above becomes O-Sea^. It is tempting to 
compare this value with that in the standard perovskite 
PbZra;Tii_j;03 that has very high efficiency for convert- 
ing electricity into mechanical strain and vice versa. Tak- 
ing for the latter an experimental "bulk" piezoelectric 
constant 633 of 11.3 C/m^ [33] and an achievable strain of 
0.2% ^24j, we obtain only 0.2eaB, four times smaller than 
the one for BN (!). Hence, the generated polarization in 
a BN sheet can be comparable or even higher than that 
in best perovskite ferroelectrics. This situation pretty 
much resembles that in the field of electrostrictive materi- 
als. Usually, the electrostriction is a small effect, because 
the induced deformation is only quadratic in electric field 
or polarization. However, dielectrics with large polariza- 
tion, such as the relaxor ferroelectric Pb(Mgi/3Nb2/3)03- 
PbTiOs, are capable of exhibiting the exceptionally large 
electrostrictive strains. 

The important characteristic, the voltage drop 
U between opposite sides of the BN strip, can be es- 
timated as (P/7reeo) ln(L/6), where P is the usual 2D 
polarization (in units of C/m), L the separation between 
the charged strip edge states (strip width), and h the 
effective radius of those states, cq the dielectric permit- 
tivity of vacuum. Taking h to be on the order of the 
lattice parameter a, for reasonable L — (100 — 1000)a, 
and P — 10~^ e/as, we formally obtain U <b — 8V for 
BN strip suspended in vacuo and smaller values for the 
strip on dielectric substrate like Si02 (e = 3.9), or more 
than 20 times that of ZnO nanowires [59]. Interestingly, 
in case of suspended BN the band bending qU exceeds the 
value of the bandgap Eg — 5.8eV [SD], meaning that it 
will result in a charge transfer between the charged edges 
to make qU « Eg. This would mean that the achieved 
bias voltage would not depend on the strip width. For 
smaller bare values of the bias, it depends on L only log- 
arithmically, i.e. very weakly. This insensitivity of the 
produced bias voltage is extremely attractive for appli- 
cations. The ability of a BN sheet to effectively convert 
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FIG. 3: Wannier functions and their shifts in a corrugated BN sheet corresponding to a zigzag period vector A(4, —4) {9 — 90°), 
A/X=0.1, phase (p— 0. (A and B) Side and top views of a vr-Iike WF contributing most to the polarization (isosurface — ±0.9). 
The yellow arrow indicate the shift of the WF in the {x, y) plane, which is antiparallel to the locally induced polarization P. (C 
and D) The (x, y)- positions of the centers of the tt and a WFs in a flat and corrugated BN sheet, respectively. The positions 
of the centers of a WFs are indicated by black squares, while those of tt by black circles. In the initial flat state, the centers of 
TT WFs coincide with those of N atoms. The red arrows show the most pronounced shifts in tt WFs induced by the corrugation 
(the amplitudes of the shifts are tripled for clarity). Note that the yellow arrow points along one of the red ones. 



wave-like deformations into electricity makes this mate- 
rial very promising for electromechanical applications at 
the nanoscale (Fig. 5). Below, we present a new con- 
cept of a nanogenerator powered by an ambient motion 
or agitation. 

The direct and efficient conversion of mechanical en- 
ergy into electricity is of great interest for practical appli- 
cations. Recently, much attention has been attracted to 
the possibility of creating miniature nano-generators us- 
ing semiconducting piezoelectric nanowires such as ZnO, 
GaN, and CdS [Ml 113 HI]; during bending they generate 
the voltage bias between the opposite sidewalls. Unfor- 
tunately, such nanowires are not really suitable for this 
application. First, as has been pointed out in Ref. [23], 
they lose a significant fraction of their generated charges 
already in the process of bending due to large internal 
conductivity by virtue of being semiconducting. Second, 
since the nanowires are grown randomly, only a small 



fraction of them (only about 1%) contributes to the gen- 
erated electric current. And third, in order to generate a 
dc current, the nanowire should slide over a metallic elec- 
trodes until the nanowire/electrode interface becomes a 
forward-biased Schottky barrier; this is equivalent to a 
ac-dc converter. 

The nanogenerator based on BN sheets does not suffer 
from all of the above shortcomings. First, the BN sheet 
is a dielectric nanostructure with a relatively wide energy 
gap (~ 5.8 eV). Second, all the atoms in the sheet are 
involved into producing electricity (it is the "bulk" ef- 
fect). Third, the nanogenerator does not need any ac-dc 
converter, like a diode bridge rectifier, because the gener- 
ated polarization has a constant direction provided that 
the direction of the wave vectors k is fixed or, in other 
words, the ambient wave-like movements are oriented in 
the space. Fourth, the produced bias voltage is nearly in- 
dependent from the BN strip width. Moreover, all these 
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FIG. 4: Wannier functions and their shifts in a corrugated BN sheet corresponding to a armchair period vector A(4, 4) {9 — 0°), 
A/A=0.1, <^= 0. (A and B) Side and top views of a a-like WF contributing most to the polarization (isosurface = ±0.9). The 
yellow arrow indicates the shift of the WF in the {x,y) plane, which is antiparallel to the locally induced polarization P. The 
{x,y)- positions of the centers of the tt and a WFs in the corrugated BN sheet (C). The positions of the centers of a WFs 
are indicated by black squares, while those of tt by black circles. The red arrows show the most pronounced shifts in a WFs 
induced by the corrugation (the amplitudes of the shifts are scaled up for clarity). Note that the yellow arrow corresponds to 
the red one. 

advantages couple with outstanding mechanical proper- 
ties of a BN sheets, like their ability to be stretched or 
bent by several percent with no atomic defects involved 

It should be stressed that the corrugation-induced po- 
larization can not be observed in a thin BN film cut out 
of the hexagonal bulk BN material. In such a film, the B 
and N arrangement is reversed in the neighboring hexag- 
onal atomic layers, so that the polarization produced by 
each such pair of layers cancels out. However, the effect 
can be observed in an ultra-thin BN films with odd num- 
ber of layers, 1, 3, 5, ... Here, of all the layers only one 
is going to produce the polarization, while the rest can 
be considered as a mechanically supporting substrate or 
scaffold. 




FIG. 5; Schematic showing generation of a bias voltage be- 
tween the edges of a corrugated BN nanosheet (the geometry 
corresponds to Fig. 4). Note that the induced polarization is 
longitudinal, because the corrugation is of the armchair type 
(see text). 
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It is worth mentioning that the BN monolayer is not 
the only 2D system capable of generating electricity un- 
der periodic bending. Its isoelectronic analog BC2N ex- 
hibits similar effect in the ground state, like, for example, 
structure found with the LDA total energy calculations 
[31]. Due to absence of inversion symmetry, BC2N in this 
structure is already polarized (electret) in the initial flat 
state with a polarization of 0.73 x lO^^e/as pointing 
along the armchair-type direction perpendicular to the 
pure C-C chains. Moreover, due to in-plane anisotropy 
under 27r/3 rotations, the corrugation-induced polariza- 
tion P (k) is no longer described by the formula ^ and 
its magnitude depends on 6 much stronger than that 
in BN. Nevertheless, there are some fc-directions along 
which the generation of polarizations is practically as ef- 
fective as in BN, such directions include the parallel and 
perpendicular ones to the above-mentioned C-C chains. 

Here, we assumed that the edges of BN sheets are 
fixed and, therefore, the surface area of the membrane 
increases due to its deformation in the third dimension, 
and this increase is simply ey. It is possible, however, 
to imagine the situation when the edges are not clamped 
and the layer bends in accordion-like fashion without any 
in-plane (average) stretching. In this case, the in-plane 
deformations dxUy and dyUx would cancel the nonlinear 
term OxUzdyU^ in Eq. (|5|, and the induced polarization 
would be defined entirely by the second, 'flexoelectric' 
term in ([s]). Such a polarization, however, is relatively 
modest because the elastic energy associated with the 
out-of-plane bending is noticeably smaller than the elas- 
tic stretching energy (this, of course, is already evi- 
dent from the comparison of the parameters a and (3). 



We have studied above the effectively infinite BN 
sheets satisfying the periodic boundary conditions, and 
the important question is whether the present results ap- 
ply to the finite 2D BN systems. This question can be 
motivated by an example of a finite graphene sheet where 
the presence of zigzag edges induces the localized elec- 
tronic states at each edge, thus modifying the properties 
of the initial infinite system j32j. Shall we expect the 
same effect in the case of BN? The answer is no, because, 
contrary to the carbon counterpart, the BN sheet does 
not have the conduction and valence bands that touch in 
a linear fashion at two inequivalcnt Dirac points, K and 
K' . This degeneracy is always lifted in a 2D BN due to a 
broken sublattice symmetry, so that the system remains 
insulating. In the present case the difference between P 
in flat and electroded corrugated states in flnite systems 
can be calculated by using periodic boundary conditions 
[33]. 

In summary, we have found that there is a strong 
polarization-corrugation coupling in a two-dimensional 
strip cut out of BN sheet, which can be decomposed 
into the nonlinear 'piezoelecric' and the flexoelectric ef- 
fects. The direction of the induced polarization strongly 
depends on the corrugation wave vector fc and changes 
approximately as (—cos 2^, sin26'), where is the angle 
between the zigzag axes and fc. The magnitude of the po- 
larization as a total dipole moment per chemical formula 
can reach very large values comparable to those in the 
best known perovskite piezoelectrics. The present effect 
may be used in nanoscale generators activated by am- 
bient vibrations that can power small electronic devices 
and circuits. 
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